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We show how one can resum the large logarithms which occur in the charm com- 
ponent of the structure function F^ c {x, Q 2 , m 2 ) (k = 2, L) taken at three light 
flavours, in the limit Q 2 m 2 . This resummation leads to an expression of the 
above quantity for four light flavours including the charm quark density. The latter 
is related to the parton densities in the three light flavour number scheme. 



The study of charm production in deep-inelastic electron-proton scattering 
has become an important issue in the extraction of parton densities in the 
proton. Thc_rjeason is that in the kinematical region covered by the HERA 
experiments ucl the charm component of the structure function i<2 )C (:r, Q 2 , m 2 ) 
constitutes about 25 % of the total structure function F%(x, Q 2 ) . Therefore 
the analysis of the data to yield parton densities can no longer treat charm 
electroproduction as a small correction. In a three flavour number scheme 
(TFNS) the charm component of the structure function can be written as 
convolutions, denoted by the symbol ®, of parton densities and heavy quark 
coefficient functions. For k = 2, L it is given by 
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where e c stands for the charge of the charm quark and the number of light 
flavours is three (nj = 3). The functions £ and A stand for the singlet 
(S) and non-singlet (NS) combinations of parton densities respectively and 
G represents the density of the gluon. Lj.^ and Hfci{i = q,g) represent the 
heavy quark coefficient functions which can also be separated into singlet and 
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non-singlet parts. The functions Lkj originate from the reactions where the 
virtual photon couples to the light quarks (u, d, s, u, d and s), whereas the Hk,i 
describe the reactions where the virtual photon is attached to the c (c) quark. 
Finally we have included in Eq. (|l|) that part of the light parton component 
of the proton structure function containing charm quark loop contributions 
to the light parton (u,d,s,g) coefficient functions only. Up to next-to-leading 
order (NLQ) the exact form of the heavy quark coefficient functions can be 
found inEl, □. If we take the limit Q 2 3> m 2 these functions behave like 
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with an analogous expression for L A ^ YMP . The asymptotic heavy quark coef- 
ficient functions have been calculated ini. When the latter are substituted in 
Eq. (|l|) we obtain the asymptotic charm component of the structure function 
which is equal to 
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A NLO analysis! reveals that for Q 2 > 20(GeV/c) 2 and x < 0.01, F 2 A C SYMP 
coincides with i?EXACT (jlj). Notice that hese values are not that far above 
threshold. However in the case of the longitudinal structure function Fl c the 
value of Q 2 is much larger and it amounts to 1000 (GeV/c) 2 for x < 0.01. 
The above findings show that in the kinematical regime covered by the HERA 
experiments the large logarithms in Eq. (|^) entirely determine Fi, c but not 
Fl,c- These large corrections vitiate the perturbation series so that they have 
to be resummed in all orders of perturbation theory. This procedure has been 
carried out in til and can be described as follows. First we add the light parton 
(u,d,s,g) component of the proton structure function denoted by Fk(3,x,Q 2 ) 
to F A C SYMP (|). Then we perform mass factorization on the asymptotic heavy 
quark coefficient functions. For H^f YMP the mass factorization can be roughly 
written as 
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where Ck,i is the light parton coefficient function appearing in -Ffc(3, x, Q 2 )- 
The quantity which absorps the charm quark mass is given by the operator 
matrix element A,-; — q,g) which is the expectation value of the heavy 
quark (here charm) operator sandwiched between light quark and gluon states. 
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Similarly one can write for Symp 

The operator matrix elements Aji. c emerge from sandwiching the light partonic 
operators between light quark and gluon states which contain charm loop con- 
tributions to the gluon self energy only. Notice that the obove relations also 
involve a redefinition of the strong coupling constant which on the righthand 
side of the above equations depencLpn four instead of on three flavours. For the 
more precise relations we refer tou. The mass factorization relations in Eqs. 
(Q) , (H) can be also used when on the lefthand side the asymptotic heavy quark 
functions are replaced by the exact ones. This has been done by the authors in 
to obtain modified coefficient functions Ck,i on the righthand side, which now 
become dependent on the charm mass. In this way one can preserve the thresh- 
old behaviour of Fk, c in their version of the variable flavour number scheme 
1VFNS) which differs from ours given below. Notice that using the method in 
there is an arbitrariness in fixing the mass dependent Ck,i- After substituting 
the mass factorization relations in the sum Fk(3,x, Q 2 ) + Fk tC (x,Q 2 ,m 2 ) one 
obtains Ffc(4, x, Q 2 ) which now depend on four flavour parton densities includ- 
ing the one representing the charm quark. The latter can be expressed into 
the original light parton densities u,d,s and g as follows 

/ c+c -(4, M 2 ) = / 4 (4, M 2 ) + / I (4,/i 2 ) 

= A% (3, 4) ® S(3, m 2 ) + A% (3, 4) ® G(3, m 2 ) , (6) 

with similar expressions for the light parton densities i. The above charm 
quark density has the property that it does not vanish at ji = m c for NLO in 
the MS scheme. Further one can construct a new charm component of the 
structure function which is expressed into convolutions of the four light quark 
densities, including charm, and the gluon density with the light parton coeffi- 
cient functions taken at four light flavours. This structure function, presented 
in the four flavour number scheme (FFNS), is given by 



pPDF 



(4, a; ,Q 2 ) = e 2 [/ c+e (4, Ai 2 ) ®C£S( 4 , 



91 

n 2 

+E( 4 , M 2 ) ®CT, s ,(4,^) +g(4, m 2 ) ®Cl g {A,z,^)] . (7) 

where C£| (purely singlet) is defined by = Cj^+C^. The above expression 
stands for the resummation of all large logarithms appearing in 

F ASYMP and 
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describes the charm structure function far-.above threshold which is equivalent 
with large Q 2 and small x. An analysis □ of the structure function F^ ACT 
(@) (TFNS) and F™ F (|) (FFNS) reveals that the former gives the best de- 
scription of charm production in the threshold region whereas far away from 
threshold it is better to use the latter structure function. Further one also 
needs a scheme which merges the advantages of these two pictures and pro- 
vides us with a good description of Fk. c in the intermediate region. This is 
given by the so called variable flavour number scheme (VFNS) for which inH 
we propose. 

F™**{x, Q 2 ,m 2 ) = F 2 P C DF (4, x, Q 2 ) + F 2 E C XACT (3, x, Q 2 , m 2 ) 

-F 2 A c SYMP (3,*,Q 2 ,m 2 ). (8) 

InJpwest order the above expression coincides with the VFNS scheme proposed 
intl. The properties of Ey™ S are discussed km and a comparison with the 
HERA data is made in From this comparison we conclude that at this 
moment the data do not allow us to discriminate between the various schemes. 
They all give a reasonable good description of the data at large as well as at 
small Q 2 . 
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